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ON  COMBINATIONS  OF  RANDOM  LOADS 

by 
D.  P.  Gaver 

and 
P.  A.  Jacobs 


Operations  Research  Department 
Naval  Postgraduate  School 
Monterey,  CA   9  3940 


1.   Introduction  and  Assumptions 

The  integrity  of  a  great  many  physical  structures  is 
potentially  threatened  by  combinations  of  physical  loads  of  vary- 
ing magnitudes  from  various  sources.   We  are  thinking  of  structures 
such  as  buildings  (for  instance  those  that  house  or  contain 
nuclear  power  plant  elements);  aircraft  and  spacecraft;  electrical 
transmission  networks;  bridges,  piers,  and  dams;  and  offshore 
oil  drilling  rigs  that  experience  loads  from  wind,  snow  and  ice, 
tides,  earthquakes,  and  so  forth.   In  many  instances  the  total 
load  or  stress  experienced  by  a  structure  varies  in  time  in  an 
apparently  random  fashion.   Certain  load  components  vary  rather 
slowly;  for  example,  that  component  resulting  from  snow  and  ice 
accumulation;  others  occur  more  nearly  as  impulses,  such 
as  those  associated  with  winds  or  earthquakes.   The  problem  is 
to  design  structures  to  withstand  the  superposition  of  loads 
from  many  sources  with  at  least  an  approximately  understood 
(high)  probability.   In  engineering  terms  we  wish  to  work  towards 
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developing  a  rational  safety  factor  technique  for  designing 
structures  to  withstand  the  combination  of  loads  anticipated. 
The  purpose  of  this  paper  is  to  describe  and  investigate  certain 
simple  but  somewhat  realistic  probabilistic  load  models  for 
use  in  design,  and  perhaps  safety,  assessment  of  structures. 

In  this  paper  we  confine  attention  to  the  superposition 
of  just  two  load  types:   shock  loads ,  and  constant  loads .   For 
example,  wind  gusts,  flash  floods,  and  earthquakes  have  varying 
magnitudes  and  have  relatively  short  durations  in  comparison  to 
the  times  between  their  occurrences;  these  will  be  modeled  as 
instantaneous  shock  loads.   On  the  other  hand  snow,  ice,  or 
water  accumulation,  or  even  the  presence  of  slow-moving  vehicles 
or  furniture,  present  loads  that  remain  nearly  constant  in  time, 
occasionally  changing  to  new  levels;  these  will  be  modeled  as 
constant  loads  that  change  infrequently.   Throughout  this  investi- 
gation it  will  be  assumed  that  the  effective  stress  exerted  by 
several  types  of  loads  acting  simultaneously  can  be  expressed 
as  a  linear  combination  (actually,  sum)  of  those  component  loads, 
the  loads  being  treated  as  stochastic  processes.   Somewhat 
similar  load  combination  models  have  been  studied  in  the  past; 
see  the  work  of  Bosshardt  [19  75],  Larrabee  and  Cornell  [19  78], 
McGuire  and  Cornell  [1974],  Pier  and  Cornell  [1973],  Wen  [1977], 
and  others.   Note  that  the  present  study  utilizes  the  notion  of 
stochastic  point  processes  as  a  modeling  tool;  see  Cox  and  Miller 
[1965]  for  an  introduction,  and  Lewis  (ed.)  [1972]  for  a  more 
extensive  treatment. 


1.1.   Modeling  Assumptions  and  Problems 

Assume  that  the  successive  magnitudes  of  shock  loads  are 
independently  and  identically  distributed  (i.i.d.)  random  variables 
with  common  continuous  distribution  function   G(y);  G(y)  =  0 
for   y  _<  0  .   The  shock  loads*  appearance  is  regulated  by  a  Poisson 
process  with  rate   p.   Furthermore,  the  constant  loads  have  i.i.d. 
magnitudes  with  distribution   F(x);  F(x)  =  0  for   x  _<  0  .   Constant 
loads  change  at  moments  of  a  Poisson  process  with  rate   A . 
Finally,  the  shock  and  constant  load  processes  are,  for  the 
most  part,  taken  to  be  statistically  independent,  although 
this  assumption  can  be  relaxed  at  times  without  difficulty. 

Let   X(t)   (respectively   Y(t))   be  the  magnitude  of  the 
constant  (shock)  load  process  at  time   t;   (Y(t)   will  usually 
be  zero).   Put   Z(t)  =  X(t)  +  Y(t),  the  superposition  of  the  two 
loads  at  time  t,  and   M(t)  =  sup   .   Z(s) ,  the  maximum  load 

combination  in   [0,t].   See  Figure  1.  Let   Ty={inf  t  >0:Z(t)  >  x}  , 

CONSTANT  LOAD  PROCESS 

x(t)T 


Y(t) 


SHOCK  LOAD  PROCESS 


Z(t) 


LOAD  COMBINATION  PROCESS 


t-*> 


Fig.  I 
THE  LOAD  COMBINATION 


i.e.    the    first-passge    time    to   a   combined   stress    level    x.    Unless    fatigue 
or   some   other  cumulative   effect   is    in   operation      T        will    represent 
the    time    to   failure   of   a    structure   whose    strength    is      x,    and   that 
is   subjected   to   a   stress   history      (Z(t),    t    >_  0}  . 

Section    2    is   concerned  with   the    description   of    the    distri- 
bution  of    the   maximum  process    and    the    distribution   of      T    .      The 
first  step   is    to   obtain  Laplace    transforms   with    respect   to    time. 
Asymptotic    results       (x   ■*■  °°)       for      E  [T    ]       and    the    distribution   of 
the   normalized  value      T*   =   T  /(E[T    ])      will   then  be    presented.       In 

XX  X 

particular  it  will  be  shown  that  the  distribution  of   T*   is 

x 

approximately   a   unit  exponential    distribution   as      x     becomes    large 
and   so   does      E  [T    ]  . 

Section    3    is    devoted   to   a    study   of   the   asymptotic 
properties   of    the    "maximum  process,"      M(t),    as      t  ■>  *>.      These 
results    relate    to    those   of   Welsch    [1972]    and   O'Brien    [1974a, b]; 
they  are    seen    to  extend   the    classical    "extreme   value"    results 
of    Gumbel    [1958]    and   Gnedenko    [1943],      We    present   results    for 
the    asymptotic    distribution   of      M(t)       and    for    the    joint   distri- 
bution  of   the    first   and   second  maximum  load    combination    in 
[0,t]      as      t  -►   oo. 


2 .  The   Laplace   Transform  of    the    Distribution   of      M(t) 

and   a    First-Passage    Time    Limit   Theorem 
2.1.       Towards    the    Distribution    of   M(t) 

Consider    first   the    distribution    function   of    the   maximum 
process,    H    (t)    =   P{M(t)     <^  x}  ,     for      x    >    0.       Its    Laplace    transform 
is   available    immediately  by  noting    that   if      T,       is    the    time    of 
the    first  change    in    the   magnitude   of    the    constant   load   process, 
then 


H     (t)     =    P{M(t)     <    x,    T       >    t)    +    P{M(t)     £   x,    T1    <    t} 


-At-       x 
e  /      exp{-yt   G(x-y)}    F(dy) 

0 

t  ,  x 

+    /      Xe  dv   /      exp{-yv   G(x-y)}    H    (t-v)    F(dy)    .        (2.1) 

0  0  x 


Next    take   Laplace    transforms   with   respect   to      t: 


CO 

h    (£)     e  /      e~Kt  H    (t)dt  =   M    (O    +    XMv(?)    h    ( O     ;  (2.2) 

X  n  X  X  XX 


reversal  of  the  order  of  integration  provides  that 


x  -       -1 

M  (O    =    I       [£  +  A  +  yG(x-y)]    F(dy)  .  (2.3) 

x       0 


Hence, 

MXU) 

V°  =  1  -  AM  (O     '  (2-4) 

x 


It  seems  to  be  difficult  to  solve  (2.1)  (invert  (2.4))  in  any 
simple  convenient  form  for  any  interesting  choice  of  the  distri- 
butions  F  and   G.   However,  useful  information  can  still  be 
gleaned  from  (2.4).   First,  it  is  clear  that  if   T    is  the 
first-passage  time  to  combined  load   x,  so 


T   =  inf{t  2  0:M(t)  >  x}  ,  (2.5) 


then,  since   P{M(t)  £  x}  =  P{T   >  t} , 


-Ct 


h  U)  =  /   e  ^   P(T   >  t)dt  ;.  (2.6) 

x      0  x 


Let      £    ■*■    0      in    (2.4)       to    find    that 


M    (0) 
m(x)     EE[Tx]    =hx(0)    =   r  .XAM    (Q)       .  (2.7) 

2\. 


We  now  record  some  expressions  for  the  mean  first-passage  time 
to  x  when  specific  distributions  for  shock  and  constant  load 
magnitudes  are  in  force. 


-x 


Example  2.1.   Identical  exponentials:   F(x)  =  G(x)  =  1  -  e   ,  x  J>  0 
In  this  case 


Mx(0)  =  x  (1  -  e  )    -   -2 


-x     -x 

ye 
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X  +  ye 
X  +  y 


-x  -i 


and 


x 


E[Tx] 


A(l-e"x)  -  yxe  x  +  ye  x  £n  [  (  X+y ),/ (  X+ye  x)  ] 


1    +   H  x   +   ^   £n[(X  +  y)/(X  +  ye    x)  ] 


(2.3) 


M-) 

y    \   x   / 


as      x   ■>   °° 


Example   2.2.       Different  exponentials 


F(x)    = 


G(x)    =   e 


-kx 


where      7-  =    2 
k 


In   this    case 


M    (0)    =    /       [X   +    ye-k(x-y)]_1  ae"ay    dy 
X  0 


=    A"2   e    2kx  N(x) 


where 


N(x)    =    A[-l    +   e2kx]    +    2y[l-ekx] 


+    2y2[-kx   +    £n[(X+y)(A    +    ye    kx)     -1]] 


Thus 

kx 


—  1  21cv  —1  e 

E[Tx]    =    A    X    N(x)      [Ae^*X   -    N(x)]    X    ~  ~-    .  (2.9) 


_    2 

Example   2.3.      F(x)    =  e    x,    G(x)    =   e         .In   this    case 


Mx(0)    =  j    [1    -   e"x]    +  H  e"x  N(x) 


where 


Thus 


N(x)    =    (e"X   -    1) [(A    +    u) (A    +    ue~X) ]_1 


E[Tx]    =   j  eX    [1    -    e    x   +    ye   x  N(x)]  [1    -    yN(x)  ]_1    -   X(\+v)    +        eX' 


(2.10) 


Example    2.4. 


F(x)    = 


1  x   <    a,  , 

-a  v      1/a 

a,  x  x  _>    a   '  , 


and 


1/a 
L2 


x   <    a"„ 


G(x) 


a2x_a  x    >    a1/* 


where        a    >    0      and      a1 ,    a2    >    0 .      In   this   case  as      x  •*•  «, 

E[Tx]    -    [Xax    +    ya2]"1x0t.  (2.11) 
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Example    2 . 5 


and 


F(x)    = 


x 


•  Ci 


G(x)     = 


x 


1 

-e 


if 


if 


if 


if 


x   <    1    , 

X     >     1     , 


X     <     1      , 
X     >     1      . 


If      a    <    3  ,    then  as      x   -> 


E  [T    ]    ~  T  xa 
x  A 


If      3    <    a,    then  as      x   -*   °° 


[T    ]    - 
x 


1     3 
—  x 

y 


(2.12) 


The  next  result  is  a  limiting  result  for  the  first 
time  the  load  combination  process  exceeds  a  given  level   x 


THEOREM  (2.1).   The  limiting  distribution  of   T    is  exponential, 
in  the  sense  that 


lim   P(m(x)  1T   >  t}  =  e  t 
x-xQ 


where 


and 


xQ  =  inf{t:F*G(t)  =  1} 


m(x)  =  E [T  ] 


Proof.      Put      p   =    vi A  and   rearrange    (2.3)    so    that   the 

denominator  of    (2.7)    is    in    the    form 

x  _! 

1    -    AM    (0)    =   F(x)    +    /      pG(x-y) [1    +    pG(x-y)]         F(dy)     =h1(x).      (2.13) 
x  0 

Now  apply  the  bounded  convergence  theorem  and  the  fact  that  F 
and  G  have  densities  to  show  that  m(x)  ->  °°  as  x  •*  xQ .  It 
next   follows  by   a   change   of   variables    in    (2.4)     that 


„  Mx(y(x) ) 

/      e    ^  P(m(x)    x   Tx    >    t}dt  =   m(x)         hx(y(x))    =   m(x)  [i-M    (Y(x))]    <2-14> 
0  x 

where   y(x)  =  £m(x)   .   By  the  bounded  convergence  theorem  applied 


to  (2.3) ,  1 

lim 

x  ->  x 


lim   Mx(y(x))  =  j      .  (2.15) 


'0 
Again  from  (2.3)  we  have 


1  -  AM  (y(x)) 


=   F(x)  +  /   [y(x)A_1  +  pG(x-y) ] [Y(x) A"1 + 1  +  pG(x-y)]"1  F(dy) 


0 


Hence , 
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m(x) [1    -    AM    (y(x)) ] 


_ 


x 


-1 


-1 


f   /       [Y(x)A         +    1    +    pG(x-y)]    x   F(dy) 
A    0 


x 


+   m(x)  [F(x)   +     /       pG(x-y)  [y(x) A    1    +    1    +    pG(x-y)  ]     X    F(dy)  ] 


Let 


x 


-1 


-1 


h9(x)    =   F(x)    +    /      pG(x-y) [y(x) A         +    1    +    pG(x-y)]    x    F(dy)    ; 
z  0 


then 


h1  (x)-h2(x) 


x 


-1 


/      F(dy)     pG(x-y) {- [1 + pG(x-y) ]  +      [Y (x) A    x + 1 + pG(x-y) ] } 
0 


x{[l    +    pG(x-y)]     [Y(x)A    1    +    1    +    pG(x-v)]}    -1. 


=    / 


x 


0 


F(dy)     pG(x-y)     y(x)A~       {[l+pG(x-y)]     [y(x)A    1  +  1+     C(x-y)]}    1 


Y(x)A    -1    k(x) 


where 


k(x)    =    / 


x 


0 


pG(x-y){[l+     pG(x-y)  ][Y(x)  A    1  +  1  +  pG  (x-y)  ]  }    1F(dy). 
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By   the  bounded  convergence    theorem 


Recall    that 


Hence , 


lim     k(x)    =0    . 
x+xQ 


j  +   0(x) 
ra(x)    =      h1(x) 


hl(x)    -   h2(x)         Ch^xJA"1      k(x) 


m 


ra~ 


A"1    +   0(x)       hl(x) 


1    +   0(x) 


k(x) 


which    tends    to   zero   as      x  ■>   xn .      Therefore, 


lim    m(x)  [1  -  AM    (y(x)  )  ] 
x  -»•  x  X 

X         XQ 


?A    -1   +    A    X     lim 


x  •>  X 


0    L 


1    + 


h2 (x)    -   h1 (x) 
h~jTx! 


U+    1)  'A    X. 


The    result  now   follows    from    (2.14),     (2.15)    and   the   unicity 
result   for   Laplace    transforms. 
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3.   Asymptotic  Results  for  the  Maximum  Load  Combination 

In  this  section  we  will  study  the  asymptotic  behavior  of 

the  maximum  load  combination  to  occur  during  time  interval   (0,t] 
as   t  -*■   °°.   Most  results  concerning  maxima  of  non-Gaussian  random 
variables  are  for  variables  in  discrete  time.   Thus  we  will  first 
obtain  results  for  an  embedded  discrete  time  load  combination 
process. 

Let   S   be  the  time  of  the  nth  change  in  magnitude  of 
the  load  combination  process;  that  is, 


S1    =  inf {t  >  0:Z(t)  t   Z(t-) } 

and  for  n  >  1  (3.1) 

Sn  =  inf{t  >  S   ,:Z(t)  ?   Z(t-)}. 


The  change  in  the  load  combination  at  time   S    may  be  due  to 

either  a  change  in  the  constant  load  process  or  to  the  arrival 

of  a  shock  load.   Let   Z   =  Z(S  )  (respectively,  X  =  X(S  ), 

n      n      r       1        n      n 

and   Yn  =  Y(Sn^  be  the  ma9nitude  of  the  load  combination 

(respectively  constant  load  and  shock  load)  at  time   S  .   Put 

n 

M«  =  0<T<n(^  •  (3-2> 

Note  that  the  times   {S.n}   are  the  arrival  times  of  a  Poisson 

process  with  rate   A  +  u   and  are  independent  of   {Z  }.   Further 

n  ' 


13 


P(ZQ  <_   x}    =    F(x) 


and  for  n  >  1 


p{Z   <_  x>  =  pF(x)  +  qF*G(x)  (3.3) 


where 


and 


00 


-At   -yt  Ja_      A 


p  =  J      Ae    e     dt  = 


0  X  +  ^ 


q  =  1-p. 


3.1.   Properties  of  An  Imbedded  Maximum  Process 

We  will  first  study  the  asymptotic  behavior  of   M   as 

n 

n  ■+■  oo.   Note  that  the  random  variables   {Z  }   are  not  independent 
In  fact   {(Xn/Zn);  n  =  0,1,...}   is  a  discrete  time  Markov  pro- 
cess.  The  following  result  describes  the  dependence  of  the 

random  variables   {Z  }. 

n 

PROPOSITION  (3.1).   The  sequence   {Z^   is  uniformly  mixing 
(cf.  Loynes  [1965,  p.  994]). 

Proof.   Let   A  and  B   be  events  such  that  A  =  f(z  ,...  Z  ) 

0  '  "       '    n 

and      B    =    g(Z         ,    Z  )  -rhon 

J     n+m'      n+m+1'    •••'•      rnen 

|P(A    n    B)     -    P(A)  P(B)  I 

=     lE[(P(B|Xn)  -    P(B));     A]  | 

=    qm|E[    /(e  (dy)     -    F(dy))    P(B|X   ^       =    y) ;    A] I 

<    2qm   P(A) 

where      ^      denotes    the   Dirac   measure    concentrated   at    the   point   x. 
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This  last  proposition  and  a  result  of  Loynes  [1965, 
Theorem   1]   imply  that  the  only  possible  nondegenerate  limiting 
distributions  of   Mn   are  the  same  three  types  that  occur  in 
the  asymptotic  behavior  of  the  maxima  of  independent  random 
variables;  namely  (except  for  scale  and  location  parameters), 


H1(x)  = 


exp{-(x~a) } 


x  £  0  , 

x  >  0,   a  >  0 


exp{-(-x«) } 


H2(x)  = 


x  <  0,   a  >  0 


x  >  0 


(3.4) 


-x 


H3(x)  =  exp{-e  A} 


<  x  < 


However,  for  the  process   {Z  },  the  following  Proposition 
(3.2)  together  with  a  result  of  Loynes  [1965,  Lemma  2] 
implies  that  the  limiting  behavior  of   M   will  not  necessarily 
be  the  same  as  that  of  the  maximum  of  a  sequence  of  independent 
random  variables  having  common  distribution   pF(x)  +  qF*G(x). 


PROPOSITION  (3.2)  .   For  any  sequences   k   ■*■  °°   and   c   ■>  °° 

2         ^  n  n 


lim   iT  X     (kn 

n  ->  °°      n  i=l 


i)  P{Z.  .  >  c  I Z1     >  c  }  >  0 
l+l     n   1     n 
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Proo 


f.   Note  that  for  any   c  >  0 


{Zi+l  >  c|Zl  >  c}  -  ql 


The  result  now  follows. 

In  order  to  obtain  the  limiting  distribution  of   Mn  we 
will  first  study  the  limiting  distribution  of  a  related  process 
which  is  defined  as  follows. 


Let 


J   =  i 

n 


1    if   Yn  =  0  , 


(   °    if   Yn 


>  0  ; 


that  is,  J   =1   if  the  nth  change  in  the  load  combination  process 
n  — 

{Z(t) ;  t  >  0}   is  due  to  a  change  in  the  constant  load  and  0  if 
it  is  due  to  the  arrival  of  a  shock  load. 


Let 


and 


T0  =  ° 


t   =  inf {k  >  t  1  :J,  =  l}.  , 
n  n-1   k 


the  index  of  the  change  in  the  load  combination  process  which  is 
due  to  the  nth  change  in  the  constant  load  process.   Let 


L  (k)  = 
n 


I 


r  .  <  n 

l  — 


tk)   i     l-l 


the   number   of      t.       that  occur  before    (discrete)     time      n      such 
that      t.    -    t._,    =   k;     (1,    | (x)    =1      if      x  =   k      and   0    otherwise) 

16 


Let 


Un  "     SUP       K    > 
xn_1£k<Tn 


the  maximum  load  combination  to  occur  during  the  nth  period  of 
time  the  constant  load  process  remains  constant. 


A 


(To) 


2  3  4  5  6 

(r.)  N       ft) 

Fig.   2 


n 


Imh  m?S^IE  time  load  combination  process 

AND  MAXIMA. 
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For  each   integer     K 

\,n    =   ^ax<n   Ui   1[0,K](Ti   "    Ti-1    }     ' 
i  — 

the  maximum  load  combination  for  the  discrete  time  process  in 

the  time  interval   [0,n]   which  is  due  to  those  shock  loads 

and  constant  loads  which  occur  during  time  intervals  of  constancy 

for  the  discrete  time  constant  load  process  which  are  of  length 

(2) 
less  than  or  equal  to   K.   Finally,  let   M     be  the  magnitude 

of  the  second  maximum  of  the  sequence 

{Ui  X[0,K](ti  "  Ti-1)?  Ti  in>  ' 

Note  from  (3.2)  that 

R,    <  M   . 
K,n  —  n 


Hence 


H    (x)  =   P(M     <  x> 
K,n  K,n  - 


K      k     Ln(k) 
=  E[  n   [F*GK(x)]  n    ] 

k=0 


>  P(M   <  x)  . 
—    n  — 


3-2-      Limiting   Results    for   Pareto-tail    Load   Distributions 
In   what    follows   we  will    assume 
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F(x)  =  axx  a  L(x)  ,      x  >  0  , 


G(x)  =  a2x  a  L(x)  ,      x  >  0 


(3.5) 


where   L   is  a  slowly  varying  function. 


Next  we  will  find  the  limitinq  distribution  of   M„     as 

K,  n 

n   ■*■    °°      for    fixed      K.      Then   we  will    show    that    the    limiting   distri- 
bution  of      M„  is    close    to    that  of      M        for      K      large. 


PROPOSITION    (3.3).       Let      F      and      G      be    as    in    (3.5),    and      u 


be    such    that 


(u    )     a    L(u    )    ~  i        as      n    ■>    -    .  (3.6) 

n  n  n 


Then 


lin  H*   .(unx)    =   exp{-x~a   P2      1      qkta!   +  ka2]}    EHk(x) 
K,n      n  ]^_q 


n->-  °° 


Proof.  Note  that 


1  -  Gk(x)  =  kG(x)  G(x)k  1 


+     (k)  G(x)2  G(x)k  2 
+  ...  +  [G(x) ]k 


ka2x  a  L(x)  as       x 
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Thus   by   result    (8.13)    on   page    278   of    Feller    [1971], 


1    -    F*Gk(x)    ~    (a1    +   ka2)x"a   L(x)     . 


Hence   as      n   -»■ 


H_.   „(u   x) 
K,n     n 


K  , 

E[exp{    I      L    (k)     in  F*GK(ux)  }] 
k=0      n  n 


K 
E[exp{    I      L    (k)  [1 
k=0 


F*G^(unx)}] 


K      L    (k) 

exp  {     I      x       [a,    +    ka„  ] 

\t0       n  X  2 


Note   that 


lim 

n   ■*■   °° 


Ln'k) 

n 


2   k 

p  q 


k   =    0,    1, . . . 


The    result   now   follows   by   the   bounded   convergence    theorem. 


THEOREM    (3.4).      Under   the   assumptions    of   Proposition    (3.3) 


-a 


lim   p{Mn  <_  unx}  =  exp{-x   [pa^^  +  qa2]} 


n  ->  co 


Proof.   Let 


K 


H  (x)    exp{  I      p2qkx  a [a,  +  ka„ ] } 
k=0  l  z 


and 
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H(x)    =     lim      H__(x) 

K  ■*■  °° 


-a 


=   exp{-x       [pa,    +   qa2 ] } 


P{M      <    u   x}    -    H(x) | 
n  —     n 


<      |P{Mn     <     UnK)     -    HK/n(UnX)l 


+     |H„       (u    x)     -    Hv(x)|     +     |H„(x)     -    H(x) 


Note    that 


|P(M^         u    x}    -    H  (u    x) |     <    E[l    -   exp{-A„(K) }]     .  (3.7) 


where 


An(K)    =    - 
n 


I        L    (k)    £n   F*Gk(u   x) 
=  K+1  n 


J         L    (k)     [1    -    F*Gk(u   x)     +   o(l    -    F*Gk(u   x)  )  ] 
k=K+l      n 


Since 


1    -    Gk(x)    <    kG(x) 


and 


1    -    F*Gk(x)     =    1  -  Gk(x)  +  k    /      G      1(y)    F(x   -   y)    G(dy) 

0 


£  k    F*G(x) 
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it    follows    that 


A    (K)    £         I         L    (k)     [k    F*G(u    x)     +    o (k    F*G(u    x))] 
n  k=K+l      n  n 


Note   that,    for      0    <    e    <    1 . 


F*G(ux)    <    F(u   x(l  -  e) )    +   G(u   x(l  -  e) )    +   F(u  xe)    G(u  xe) 
n        —  n  n  n  n 


Therefore  for    fixed     6     >    0      there    is    an      Nfi      such    that    for 
n    >  NQ 

00  °°      L    (k) 

1  L    (k)    k   F*G(u   x)    <         I      -^ —  k[a   x        +   a?x        +  6  ] 
k=K+l  n                           n             k=K+l      n                x                  z 


Further, 


n 


2    k 


E[Ln(k)]    <    E[    I      l{k}(T.    -    T._1)]    =    npV    . 

i=l 


Hence    for  n    >_  Nfi(5  ) 


-ot    „  ,     r      2  k, 


E[  I  L  (k)  k  F*G(u  x)]  <_    [a,  +a9]x  a  +  6  ]    J    pqk.   (3.8 

k=K+l    n  n  i         Z  k=K+i 

Choose  y  >  0.   Apply  Jensen's  inequality  to  (3.7): 


for   n  >  NQ  (5  ) 


E[l    -    exp{-A    (K)}]    <    l-exp{E[-A    (K)  ]  }    ■*    0 
n  —  c  n 
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as      K   +   °°.      Thus    for   sufficiently   large      n      and     K 


|p{Mniunx}    -  HK>n(unx)l 


is    arbitrarily   small.       It   follows    that 


lim     P{Mn    <_  xun)    =   exp{-x~°t[pa1    +   qa2]}  (3.9) 


whereas,    if   the   load  process  were   one   of   independent  random  vari- 
ables,   denoted  by      Z'      with    distribution      pF(x)    +   qF*G(x),    the 
maximum  would  be    distributed  as    follows 


lim      P{M'     <    xu    }    =   exp{-x        a,    +   qa„ ]}     ;  3.10 

n   —        n  *  1        ^    2 

n   -*■   °° 


The      u        is    the    same    in   each   case;    see    (3.6).      Comparison   of 
(3.9)    and    (3.10)    shows    that   the  expression    (3.10)    overestimates 
the    probability   of  exceeding   a   given   stress    level. 

O'Brien    [19  74a]    obtained   the   above   result   in    the   case 
in  which   there    is   a   constant   load    (i.e.,    a2   =   0)      alone,    or 


only   a   shock    load    (i.e.,    a,    =    0)  . 
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3.3.      Results    for   Continuous   Time 

THEOREM    (3.5).      Let      F      and      G      be    as    in    (3.5).       Let      u(t) 
be    such    that      u(t)~a   L(u(t))~  1/t      as      t  •*   <*>.      Then       A+y. 


lim      P(M(t)    <    xu(t)}    =   exp{-x~a[Aa,    +    ya„]}    .         (3.9) 


Proof.       Let      N(t)       be    the  number  of   changes    in   the   load   combina- 
tion  process,    (Z(u);    u    >_  0},    in    the    time    interval       [0,t).      The 
process      {N(t);t  _>    0}      is    a   Poisson   process   with    rate 
Further, 

M(t)    =   H,(t)     . 

The  result  now  follows  from  the  proofs  of  Proposition  (3.3) 
and  Theorem  (3.4)  ,.  the  strong  law  of  large  numbers  for  N(t), 
and  the  bounded  convergence  theorem. 

Slight  modifications  of  the  proofs  of  Proposition  (3.3) 


and  Theorems  (3.4)  and  (3.5)  imply  the  following  result. 


THEOREM.  (3.6)  ,   Let 


F(x)  =  x"aL(x)    and    G(x)  =  x~6L(x) 


where   L   is  a  slowly  varying  function. 


a)   If   a  <  6   and   u(t)   is  such  that   u(t)  aL[u(t) ]  ~  i/t   as  t 

then 

lim  P(M(t)  <_  xu(t)>  =  exp(-Ax   }.  (3.10) 

t  -*■  °° 
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b)   If  3  <  a   and   u(t)   is  such  that   u(t)  6L[u(t)]  -  1/t  as   t  +   <* 

then 

lim  p(M(t)  <_  xu(t)  }  =  exp{-  yx  P  }  . 

t  ■*■  °° 

The  above  result  indicates  that  if  the  tail  of  the  dis- 
tribution of  the  magnitude  of  a  constant  load  (respectively  shock 
load)  dominates  that  of  the  shock  (respectively  constant  load) , 
the  asymptotic  behavior  of  the  maximum  load  combination  is  the 
same  as  that  of  the  maximum  constant  load  (respectively  shock 
load)  by  itself.   On  the  other  hand,  Theorem  (3.5)  indicates 
that  if  tails  of  the  distributions  of  the  magnitudes  of  the 
constant  and  shock  loads  are  comparable,  then  the  asymptotic 
behavior  of  the  maxima  of  the  load  combination  process  depends 
on  both  the  constant  and  shock  load  processes. 

3.4 .   Limiting  Results  for  the  Joint  Distribution  of  the  First 

and  Second  Maxima 

(2) 
Let   M      denote  the  magnitude  of  the  second  maxima  of 
n 

{Z,  ;  k  <_  n}.   Assume   F   and   G   are  as  in  (3.5).   Let 


H(x)  =  exp{-x   [pa,  +  qa~  ]  } 


as  before.   Techniques  similar  to  those  used  in  the  proofs  of 
the  results  in  Subsection  (3.2)  yield  the  following  result 
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THEOREM  (3.7) .   Let   u   be  as  in  (3.6) .   Then 

(2) 


lim  PtM;    -  yun'  l\i   -  xun}  =  B(^'x) 


n  ■*■   °° 


where  for   x  <  y, 


B(y,x)  =  H(x) 


and  for   y  <  x 


B(y,x)  =  H(y)  {1  +  [p ^  +  qa2J[y  a  -  x~a]}  . 

A  result  of  Welsch  [1972]  implies  that  the  limiting 
distribution   B   must  be  of  the  following  form 


B(y,x) 


H(x) 


if  y  >  x, 


(3.1] 


H(y){l-g[(£n  H(x)/£n  H(y))]£n  H (y) }  ,   if   y  <  x 


where   g(s),  0  <  s  <_  1   is  a  concave,  nonincreasing  function  which 
satisfies   g(0) (1  -  s)  <  g(s)  <  1  -  s.   In  our  case 


g(s)  = 


p  al  +  qa2 
pax  +  qa2 


(1  -  s)  . 


(3.12) 
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(2) 
Let      M         (t)       denote    the   magnitude   of    the    second  maximum 

to   occur   in       [0,t]       for    the    continuous    time    load   combination 

process      {Z(s);    s   ^  0}.      Arguments    similar    to    those    in   Section 

3.3   can  be    used   to   obtain    the    following   result. 


THEOREM    (3.8).      Let      F      and      G     be    as    in    (3.5)    and      u(t)      be 


-a 


such    that      u(t)         L(u(t))    ~   1/t      as      t 


For      x    >    y 


lim   P(M(t)     <    xu(t),    M(2)(t)     <yu(t)} 

t     ->     00 


=  exp{-y       [Aa1    +    ya2] } 


-a 


■a. 


1  +  A+ii  ai +  ya2  {y    -  x  ) 


for      x    <    y 


(2) 


lim      P{M(t)    <    xu(t) ,    Mv    ; (t)     <    yu(t)} 


t   -»• 


=   exp{-x       [Xa,    +      ya?]}. 
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